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The algebraic formulation of Wick’s theorem that allows one to present the vacuum or thermal 
averages of the chronological product of an arbitrary number of field operators as a determinant 
(permanent) of the matrix is proposed. Each element of the matrix is the average of the chronological 
product of only two operators. This formulation is extremely convenient for practical calculations 
in quantum field theory and statistical physics by the methods of symbolic mathematics using 
computers. 

PACS numbers: 02.10.Ud, 02.70.Wz, 05.30.Fk, ll.lO.Wx 


INTRODUCTION 

Wick’s theorems are used extensively in quantum field 
theory and statistical physics ■ They allow one 
to use the Green’s functions method, and consequent 
to apply the Feynman’s diagrams for investigations [ll-lsi . 
The hrst of these, which can be called Wick’s Theorem 
for Ordinary Products, gives us the opportunity to reduce 
in almost automatic mode the usual product of operators 
into a unique sum of normal products multiplied by c- 
numbers. It can be formulated as follows [^. Let Ai {xi) 
{i = ) are “linear operators”, i.e., some lin¬ 

ear combinations of creation and annihilation operators. 
Then the ordinary product of linear operators is equal to 
the sum of all the corresponding normal products with all 
possible contractions, including the normal product with¬ 
out contractions, i.e., 

Ai... An = : Ai... An ■ '. A 1 A 2 ... An: -I-... 

I_I 

-b: Ai ... An-i An : -b: Ai ... An : 

I_I I-1 

-b : A 12 I 2 bl3^4 ... : -b •.. , 

I_II_I 

where AiAj = AiAj— : AiAj : (*, j = 1, 2,..., n) is the 

I_,_I 

contraction between the factors Ai and Aj. Since the 
vacuum expectation value of the normal ordered product 
is equal zero, this theorem provides us a way of expressing 
the vacuum expectation values of n linear operators in 
terms of vacuum expectation values of two operators. 

Wick’s Theorem for Chronological Products asserts 
that the T-product of a system of n linear operators is 
equal to the sum of their normal products with all possi¬ 
ble chronological contractions, including the term without 
contractions. It follows directly from the previous theo¬ 
rem and gives the opportunity to calculate the vacuum 
expectation values of the chronological products of linear 
operators. 

Finally, from Wick’s theorem for chronological prod¬ 
ucts the Generalized Wick’s Theorem can be ob¬ 
tained. It asserts that the vacuum expectation value 
of the chronological product of n -\- 1 linear operators 


A,Bi,..., Bn can be decomposed into the sum of n vac¬ 
uum expectation values of the same chronological prod¬ 
ucts with all possible contractions of one of these opera¬ 
tors (for example A ) with all others, i.e., 


{TiAB,...Bn))o= X {t (aB, ... Bi .. Bn)) . 

l<i<n ° 

,_, ( 1 ) 

Here AA = T (AiAj) - : A^Aj := {T{AiAj))^ {i,j = 
1,2,... ,n ) is the chronological contraction between the 
factors Ai and Aj. It should be noted that, in contrast 
to the usual Wick’s theorem for chronological products, 
there are no expressions involving a number of contrac¬ 
tions greater than one on the right-hand side of ©• 
Wick’s theorem for chronological products or its ge¬ 
neralized version are used for the calculation of matrix 
elements of the scattering matrix in each order of pertur¬ 
bation theory The procedure is reduced to calcula¬ 

tion of the vacuum expectation of chronological products 
of the field operators in the interaction representation. 
As factors in these products a number of operators ifi 
of the Fermi fields and the same number of their con¬ 
jugate operators 'ifi , as well as operators of the Bose 
fields ips = ^ may be used. Here all continu¬ 

ous and discrete variables are included in the index. In 
the interaction representation the operators and 

correspond to free helds and satisfy the commuta- 


tion relationships of the form 

i’i, il’j] 

+ = [V’i, ^0 

['0i ) V’j ] ^ — ^ij 5 

pi \ pi ^ 

= 





= Drs- Therefore, the averaging of the 


Fermi and Bose fields can be performed independently. 


BASIC THEOREM 

Since we may rearrange the order of the operators in¬ 
side T-products taking into account the change of the 
sign, which arises when the order of the Fermi operators 
is changed, we present our vacuum expectation value of 
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the chronological product of the Fermi operators in the 
form 

± {T (V’»2^i2) • ■ • • (2) 

To calculate (I2|), we can use Wick’s theorem for chrono¬ 
logical products. However, while considering the higher- 
order perturbation theory, the number of pairs ipiijjj of 
operators ipi and ipj becomes so large that the direct 
application of this theorem begins to represent certain 
problems because it is very difficult to sort through all 
the possible contractions between 'ipt and ipj. 

A consistent use of generalized Wick’s theorem would 
introduce a greater accuracy in our actions. However, in 
this case we expect very cumbersome and tedious calcu¬ 


lations. Hereinafter we show that the computation of ([2]) 
can be easily performed using a simple formula 

(T [{t/ji.ijj,) ... )o = det J , (3) 

where 

I-^—I 

(a,^ = l,2,...,n) . (4) 

The proof of this theorem is by induction. Let us assume 
now that is true for n pairs ; and consider it 

for the case n -|- 1. Using generalized Wick’s theorem we 
have 


{T {^ 2 ^ 02 ) ■ ■ ■ + [(V’in.lV’il) (V’i 2 ^^ 2 ) ' ' ' 

n—1 

7=2 

^in + lJn + l 


Taking into account ® we obtain 
J 


(T [{iJi.tljj,) = 


^hjn + l 

^i2jl 

^*n + li2 
^^232 

^in + ljn. 
^i23n 

^i23n + l 


^ilj2 
^in + lj2 

■ ■ 

^in + ljn 


^injl 

^in,32 

■ ^injn 


^injl 

^inh 

■ ■ 


^ii3i 

^ilt2 

■ ^hjn 


^iljl 

^ilj2 

■ 




+ ^in + ljn + 1 

^i2jl 

^i232 

^*2jn 


^in-lj2 






^in + lt2 

^in + l3n 


^injl 


■ ^injn 


Rearranging the rows in the determinants in ([5]) it is easy 
to see that the right hand side is the expansion of the 




^iljl 

^il32 

^UJn + 1 



^*2^1 

^1232 

^i23n. + l 

det 

— 






^in + ljl 

^in + l32 ■ ■ ■ 



(a,/3 = 1, 

2, . . . , 7Z + 1) 



along the last column Q. The validity of (|31) for n = 1 
follows from the definition in Q. 


Note that this result does not depend on the way 
how we divide the operators on the left hand side of ([3]) 
into pairs Indeed, if on the left side of ([3]) we 

permute, for example, and ipj^ (/3 ^ 7 ), it changes 
its sign. The same happens on the right hand side of 
since this change leads to the permutation of two 
columns in the determinant, and it also changes its 
sign. Similarly, in the case of a permutation of and 
'i’is Obviously, when the whole pair ipia4’jf3 is 

transposed, the left and right hand sides of ® do not 
change. 
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Obviously, the formula similar to ( 3 ) can be obtained 
and in the case of Bose fields 


(t [( 4 ^ Vi (ViVi 

= perm(Ai^j,3) , 


where 


= ViVi ^ = (t (ViVi 

(a,/ 3 = l, 2 ,...,n) . 


( 6 ) 


( 7 ) 


In quantum statistics the n-body thermal, or 
imaginary-time, Green’s functions in the Grand Canoni¬ 
cal Ensemble are defined as the thermal trace of a time- 
ordered product of the fields operators in the imaginary- 
time Heisenberg representation [m. To calculate them 
in each order of perturbation theory, Wick’s theorem is 
also used. Obviously, in this case, the theorem also may 
be formulated in the form (jS)) and (jG]) convenient for prac¬ 
tical calculation. 


where V (xi — X2) = U {vi — r2) 5 (ti — T2). We can im¬ 
mediately take into account the reduction of the dis¬ 
connected diagrams, if we assume in ([8]) {S)q = 1 and 


A (xi, 0:2) = 0 [6|. Then 


1/7 (xi,X2) =-A (xi,X2)-I-2^^ J dzidz 2 V {zi - Z 2 ) 
X A(xi,Z2 ) 


-A(xi,zi) 


A(zi,X2) A(zi,zi) 
A(z2,X2) A{z2,Zi) 
A(zi,X2) A{zi,Z 2 ) 
A{z2,X2) A(z2,Z2) 
= -A (xi,X2) 


-I- j dzidz2A (xi, 0i) V {zi - Z2) A (^i, ^2) A (2:2, X2) 
- J dzidz2A{xi,zi)V {zi -22) A (22,-22) A (21, X2) . 


SIMPLE EXAMPLES 


Taking into account (xi,X2) = —A(xi,X2) we ob¬ 
tain finally 


In order to demonstrate the usability of the proposed 
formulation of Wick’s theorem, we find the first-order 
correction to the one- and two-particle thermal Greens 
functions for the Fermi system described in the interac¬ 
tion representation by the Hamiltonian 

H^nt ('^) = ^ J dridr2tpa (ri,r)V (r2,r) 

xU (ri -r2)ipi3 (r 2 ,r) ipa (ri,T) 

that contains the product of the field operators ipa (ri, r) 
and ipa (ri,T) in this representation (parameter a indi¬ 
cates the spin projections, r is the imaginary-time). The 
one-particle Green’s function can be represented as 

(TV V {xi)tpix2)S])^ 


Ql (xi,X2) = -- 


(>5)7 


where (.. .)q is the symbol for the Gibbs average over 
the states of a system of noninteracting particles, x = 
(r,T, a) and 

S {t) = Tr exp I - y dr' Hint [t') 

We obtain 

5 /(xi,X 2 ) = -{Tr V {,Xi)'lp{x 2 )])^ 

+ J dzidz2V (21 - 22) 

x{Tr [{ip {xi)Pp (X2)) {ip (21) Ip (21)) {ip (22) Ip (22))] )q 
= -A (xi, X2) -I- 2(^ j ^^^d,Z2V (21 - 22) 


A(xi,X2) A(xi,2i) A(xi,22) 
A(2i,X2) a (21,21) A (21,22) 
A(22,X2) a (22,21) A (22,22) 


( 8 ) 


Gl (xi,X 2 ) = G^°^ (xi,X 2 ) 

-b J dzidZ 2 G^°Hxi,Zi)i:'^~^'> {zi,Z 2 ) 2 G^°Hz 2 ,X 2 ) , 


(2:1, -22) = -V (21 - 22) (21,22) 

-1-5(21 — 22) J dzV {zi — z) G^°^ (-2,2) . 


Similarly, for the two-particle Green’s function 


Gn (X 1 ,X 2 ,X 3 ,X 4 ) = 


(Tr [ip {xi)ip (xa) Ip (X3) Ip (X4) S] 

( 5 )o 


we have 
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2 (‘ 5 )c 


Qll {xi,X2,X3,Xi) = (T [{-Ip (xi)-Ip {xs)) {'Ip{x 2 )'lp{x 4 ))])g 
dzidz2V (zi - Z2) {T [{pj (xi) ^ (xa)) {ip {x2)'ip (0:^4)) {ip (21) V' (^i)) {tp (2:2) V' (22))] )(- 


A(xi,a;3) A(xi,a:4) 
A(x2,a;3) A(a:2,a:4) 


0 0 A(a;i,zi) A(a;i,Z2) 

0 0 A (3:2,21) A(x 2 ,Z 2 ) 

A (01,3:3) A (01,3:4) A (01,01) A (01,02) 

A (02,03) A (02,3:4) A (02,01) A (02,02) 

= (cCi, X3) (3:2, 3:4) - {xi,X 4 ) G^°^ {x2,X3) 


j dzidz2V (01 — 02) 


- / dzidz2G‘'^'’ (a:i,0i)^^°^ (3^2,-22) V (01 - 02) (21,0:3)0^°^ (22,0:4) - (21,3:4)^*^°^ (22,^3) 


CONCLUSIONS 


The above theorem has an important consequence. In 
fact, it establishes a perfect coincidence between the vac¬ 
uum expectation values of the chronological products of 
n pairs of field operators and the n-order determinant. 
If to present this determinant as the sum of the elements 
and cofactors of one any row or column, and thereafter to 
use again the indicated coincidence for the (n — l)-order 
determinants included in each summand, we will return 
to the generalized Wick’s theorem. Alternatively, we can 
select in the our n-order determinant arbitrary m rows or 
columns (I < to < n) and use the Generalized Laplace’s 
Expansion Q for its presentation as the sum of the prod¬ 
ucts of all TO-rowed minors using these rows (or columns) 
and their algebraic complements. Then, taking into ac¬ 
count our theorem, we obtain a representation of the vac¬ 
uum expectation values of the chronological products of 
n pairs of field operators as the sum of the products of 
vacuum expectation values of the chronological products 
of TO pairs of operators and vacuum expectation values 
of the chronological products of n — to pairs. The num¬ 
ber of terms in this sum is equal to n\/ml (n — to)!. This 
decomposition can be useful for the summation of blocks 
of diagrams. 


Representations ([ 21 ) and (O not only greatly simplify 
all calculations, but also allow one to perform them using 
a computer with programs of symbolic mathematics Q. 
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